Anisotropic models of dark energy dominated Universe 



43 



Of 



Suresh Kumar * 
July 26, 2012 

Abstract 

In this paper, we investigate the role of dynamical dark energy in the evolution of Universe within 
the framework of a spatially homogeneous Bianchi type-I space-time filled with non-interacting matter 
and dark energy fluids. We utilize a special law of variation for the average Hubble's parameter in 
Bianchi-I space-time that leads to power-law and exponential-law cosmologies. Also, it yields constant 
value of deceleration parameter given by q = n — 1, which provides accelerating models of the Universe 
for n < 1 and decelerating ones for n > 1. In the present work, we emphasize that decelerating models 
can be described by the usual matter fluid, while dynamics of accelerating Universe can be described by 
(3JT)[ considering some exotic type of matter such as the dark energy. We study the role of dark energy with 

variable equation of state parameter in the evolution of power-law and exponential-law cosmological 
models. The analysis of the models reveals that the present acceleration, isotropy and flatness of the 
Universe are the natural consequences of dark energy. The equation of state parameter of dark energy 
is found to be consistent with the seven year WMAP observations. The analysis of the lookback time, 
43 . proper distance, luminosity distance, angular diameter distance and event horizon is also carried out 

for the power-law and exponential-law models. 

Keywords: Bianchi-I space-time • Average Hubble's parameter • Accelerating Universe • Dark energy 
PACS number: 98.80.-k 



1 Introduction 

The late time acceleration of the Universe is confirmed by different data sets of complimentary nature 
such as type la supernovae (SN la), baryon oscillations, galaxy clustering, cosmic microwave background 
(CMB) and weak lensing [TJ. The cause of accelerated expansion of the Universe has been attributed 
to some unknown form of energy named as the dark energy (DE). It is a hypothetical form of energy 
that permeates all of space and tends to increase the rate of expansion of the Universe [2|. The most 
recent WMAP observations indicate that DE accounts for 72% of the total mass energy of the Universe 
|3j. Many cosmologists believe that the simplest candidate for the DE is the cosmological constant (A) 
or vacuum energy since it fits the observational data well. During the cosmological evolution, the A-term 
has the constant energy density and pressure p( de ^ = —p( de \ where the superscript (de) stands for DE. 
However, one has the reason to dislike the cosmological constant since it always suffers from the theoretical 
problems such as the "fine-tuning" and "cosmic coincidence" puzzles [3]. That is why, the different forms 
of dynamically changing DE with an effective equation of state (EoS), w^ de ' = p( de )/p( de ^ < —1/3, have 
been proposed in the literature. Other possible forms of DE include quintessence (w^ de ^ > — 1)[5], phantom 
( w (de) ^ _-q ^gj e ^- c rp^g observational constraints on the EoS parameter have been found in various 
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works in literature. For instance, the seven year WMAP results put w^e in the range —1.55 < Wd e < —0.7 
(see, Jarosik et al. [7]). 

Even though the Universe on a large scale, appears isotropic and homogeneous at the present time, 
there are no observational data that guarantee in an epoch prior to the recombination. The sorts of 
matter fields in the early Universe are uncertain. In the early stages of evolution, the Universe could not 
have had such a smoothed out picture because near the big bang singularity neither the assumption of 
spherical symmetry nor of isotropy can be strictly valid. Therefore, to describe the evolution of Universe 
right from the big bang, it appears appropriate to suppose a geometry that is more general than just the 
isotropic and homogeneous FRW geometry. In this regard, Bianchi type-I space-time is of fundamental 
importance since it provides the requisite framework. A spatially homogeneous and anisotropic Bianchi 
type-I model is considered as the simplest generalization of the FRW flat model. It is described by the 
line element (in units c = 1) 

ds 2 = -dt 2 + A 2 dx 2 + B 2 dy 2 + C 2 dz 2 , (1) 

where A , B and C are the metric functions or directional scale factors of cosmic time t. If any two of 
the directional scale factors are equal and third one is different (say, A ^ B = C), the space-time is said 
to be axially symmetric or locally rotationally symmetric (LRS). In case A ^ B ^ C, the space-time is 
totally anisotropic. 

Following Berman [8], Kumar and Singh [9 J proposed a special law of variation for the average Hubble's 
parameter in Bianchi-I space-time, which yields a constant value of deceleration parameter (DP). The 
law reads as 

H = D(ABC)-? = Da~ n , (2) 
where D > 0, n > are constants and a = (ABC)s is the average scale factor of the Bianchi-I space-time 

Such a law of variation for Hubble's parameter is not inconsistent with the observations and is also 
approximately valid for slowly time- varying DP models. The law provides explicit forms of scale factors 
governing the Bianchi-I Universe and facilitates to describe accelerating as well as decelerating modes of 
evolution of the Universe. Models with constant DP have been extensively studied in the literature in 
different contexts (see, Kumar and Singh [9] and references therein). Most of the models with constant 
DP have been studied by considering perfect fluid or ordinary matter in the Universe. But the ordinary 
matter is not enough to describe the dynamics of an accelerating Universe as mentioned earlier. This 
motivates the researchers to consider the models of the Universe filled with some exotic type of matter 
such as the DE along with the usual matter fluid. Recently Akarsu and Kilinc |1Q[ [TTT 112] . Yadav et al. 
P31 El E5], Pradhan et al. [HI E], Adhav et al. [HIS EH] and Kumar et al. [2H EH E3] have studied 
the role of DE with variable EoS parameter in cosmological models evolving with constant DP. 

In this paper, we have considered non-interacting matter fluid and DE energy components with con- 
stant DP within the framework of the Bianchi-I space-time JI]). The EoS parameter of DE has been 
allowed to vary with time. The paper is organized as follows. In Sect. 2, the field equations have been 
presented. Sect. 3 deals with the solutions of the field equations and physical behavior of the model. Sect. 
4 is devoted to the study of some observational quantities such as lookback time, proper distance, lumi- 
nosity distance, angular diameter distance and event horizon for the derived models. Finally, concluding 
remarks have been given in Sect. 5. 

2 Field equations 

The Einstein's field equations in case of a mixture of perfect fluid and DE components, in the units 
8ttG = c = 1, read as 

Rij ~ ~^9ijR = (3) 
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where Tjj = T^J 1 + T^ 1 is the overall energy momentum tensor with Z^- and T^- as the energy 
momentum tensors of matter and DE, respectively. These are given by 

T { f { = diag [-p {m \ p {m) ,p (m) , p (m) ] 

= diag [-1, w {m) ,w {m \ w {m) ]p {m) (4) 



and 



T {de)i = diag [-pV*\ J de \ pW] 



diag [-1, ,w ( - de \ w^p^ (5) 



where p( m ) and p^ are, respectively the energy density and pressure of the matter fluid component while 
w (m) — p(m) J p{m) - g ^ s -g Q g p arame ter. Similarly, p( de ^ and p( de ^ are, respectively the energy density and 
pressure of the DE component while u/ de ) = j)( de ) / p( de ) is the corresponding EoS parameter. 

In a comoving coordinate system, the field equations ([3]), for the Bianchi-I space-time (pQ), in case of 
flU) and ©, read as 

f + 2. + = - w ( m )p( m ) _ „,(*)„(*), (6) 

(_y .ri (_y .rl 

d + | + 4| = - a ,<»V»> - „,<*),,<*), (8) 

Here the over dot denotes derivative with respect to t. 
The energy conservation equation T l3 ,- = yields 

p(m) + 3(1 + w (™)} p (m) H + ^(de) + 3^ + w (de)) p (de) H = ^ (10) 

where -ff = d/a. 

3 Solution of Field Equations 

Subtracting ([6]) from ([7]), ([6]) from ([8]), ([7]) from ([8]) and taking second integral of each, we get the following 
three relations respectively: 

— = di exp ^xi J a~ 3 dt^J , (11) 
^ = d 2 exp ^x 2 J a~ 3 dt \ , (12) 
^ = d 3 exp (2:3 / a' 3 dt] , (13) 



C 

where d\ , x\ , d 2 , x 2 , d^ and X3 are constants of integration 
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From (jlip -(|13 p . the metric functions can be explicitly written as 



A(t) = aiaexp (bi J a~ 3 dt\ , (14) 

B(t) = a 2 aexp (b 2 J a' 3 dt \ , (15) 
C(t) = a 3 aexp J a~ 3 dt \ , (16) 



where 



ai = y / d 1 d 2 ,a 2 = a d 1 1 d 3 ,a 3 = ^/ (d 2 d 3 ) r , 



Xl+X 2 , X3-X1 -(X 2 +X 3 ) 

61 = ^— ,6 2 = ^— ,63 = § 



These constants satisfy the following two relations: 

aia 2 a 3 = 1, bi + 6 2 + 63 = 0. (17) 

The field equations ©-dH]) involve seven unknown variables, viz., A, B, C, w^ rn \ w^ de \ p( m ) and 
p( de ) . Therefore, to find a deterministic solution of the equations with the law ([2]), we need two suitable 
assumptions connecting the unknown variables. 

Following Akarsu and Kilinc [10], we assume that the matter fluid and DE components interact 
minimally. Therefore, the energy momentum tensors of the two sources may be conserved separately. 

The energy conservation equation % \- = 0, of the matter fluid leads to 

pW+3(l + t»W)pWff = 0. (18) 

Therefore, (ITU]) gives 

the energy conservation equation of DE component. 

Next, we assume that the EoS parameter of the matter fluid to be a constant, that is, 

w = = const., (20) 

while has been allowed to vary with time since the current cosmological data from SN la, CMB 

and large scale structures mildly favor dynamically evolving DE crossing the phantom divide line (PDL) 
( w {de) = -1) [TUJ. 

In view of ([2U]) . integration of ([T8]) leads to 

pM =Coa - 3 ( 1 +- (m) ), (21) 

where cq is a positive constant of integration. 

In the following subsections, we discuss the DE cosmology for and n = by using the law ([2]). 
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3.1 DE Cosmology for n^O 

In this case, integration of ([2]) leads to 



where c\ is a constant of integration. 
Inserting (|22j) into JH])-JT6]), we have 



a(t) = (nDt + ci)' 



(22) 



A(t) = a\ (nDt + c\) « exp 

Bit) = a2(nDt + ci)« exp 
C(t) = a^(nDt + ci) « exp 



6i 


D(n- 


co 
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D(n- 


CO 






D(n- 


CO 



(nDt + ci)' 



(nDt + ci) 



(nDi + a) 



(23) 

(24) 
(25) 



provided 

Therefore, the model (pQ) becomes 



(26) 



where T = nDt + c\ and = D ^s) ■ 

The Hubble parameter (H), energy density (p^ 1 ) of perfect fluid, DE density (p( de ^) and EoS param- 
eter (ti/ de )) of DE, for the model (J2SD are found to be 



p (m) =co(n,Dt + ci)" 



i? = DfnDi + ci) -1 , 

■ 3(l+iu( m )) 



-3(lW">) 



i/7 



(de) 



P 



(de) 



(2n - 3)D' z (nDt + a)' 2 - p {nDt + a)^ - c w (m> [nDt + a 



-3(lW m >) 



(27) 

(28) 
(29) 
(30) 



where j3 = bf + bl + &i& 2 - 

The above solutions satisfy the equation (fT9|) identically, as expected. 
The spatial volume (V") and expansion scalar (9) of the model read as 

V = a 3 = (nDi + ci)™, 



= 3H = 3(nDt + ci)" 1 . 
The anisotropy parameter (.A) and shear scalar (cr) of the model are given by 



A 



1 



9H 2 



A-b\ 2 (c_aX 

a~b) + \b~cI + \c~a) 



^(nDi + ci) n 



(31) 
(32) 

(33) 
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a 2 = ^AH 2 = PoinDt + C!)^. 
The density parameter of perfect fluid and the density parameter Q,^ de ^ of DE are as follows: 

->( m ) fn 2n-3(l + m < m )) 



CQ 



3# 2 3D 2 



(nDt + ci) 



(34) 



(35) 



n (cfe) 



1 



3tf 2 D 2 
The value of DP (g) is found to be 



(nDt + ci) 



-2(3-n) 



CQ 

3D 2 



(nDt + ci) 



2n-3(l+m< m )) 



aa 

T2 = "-1, 
or 



(36) 



(37) 



which is a constant. A positive sign of q, i.e., n > 1 corresponds to the standard decelerating model 
whereas the negative sign of q, i.e., < n < 1 indicates acceleration. The expansion of the Universe at 
a constant rate corresponds to n = 1, i.e., q = 0. Also, recent observations of SN la (|24J-|31J) reveal 
that the present Universe is accelerating and value of DP lies somewhere in the range — 1 < q < 0. In a 
recent paper, Kumar [32] found DP in the same range for the power-law cosmology by using latest data 
from SN la and H(z) observations. It follows that in the derived model, one can choose the values of DP 
consistent with the observations. 




Figure 1: Scale factors vs time with b\ 
as unity. 



0.3, b 2 = 3.3, b 3 



-3.6, D = 2, n = 0.5 and other constants taken 



We observe that at t = —c\/nD, the spatial volume vanishes while all other parameters diverge. 
Therefore, the model has a big bang singularity at t = —ci/nD, which can be shifted to t = by choosing 
ci = 0. The singularity is point type as the directional scale factors A(t), B(t) and C(t) vanish at the 
initial moment. The cosmological evolution in Bianchi-I space-time is expansionary, since A(t), B(t) and 
C(t) monotonically increase with the cosmic time as may be observed from Fig.l. Also, we see that the 
scale factors start growing at different rates in the beginning. This confirms the anisotropy of the derived 
model in the early Universe. But later on, the growth of scale factors takes place at the same rate. So the 
Universe expands at same rate in the three spatial directions, i.e., the Universe becomes isotropic. Since 
the present-day Universe is isotropic, so the derived model is consistent with the observations. 
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t 

Figure 2: «;(*) vs t with n = 0.05, D = 2, /3 = c = ci = 1, w {m) = 0. 

The parameters H, p^ m \ p^ de \ 6 and a 2 start off with extremely large values, and continue to decrease 
with the expansion of the Universe. The anisotropy parameter A also decreases with the cosmic evolution 
provided n < 3. This shows that anisotropy of the model goes off during the cosmic evolution. The spatial 
volume grows with the cosmic time. The EoS parameter of DE asymptotically approaches to — 1 
provided re < 3. Fig. 2 depicts the variation of u/ de ) during the cosmic evolution for some selected values 
of constants. We observe that w^ de ^ starts in phantom region {w^^ < — 1), crosses the PDL (ti/ rfe ) = —1) 
and finally varies in the quintessence region 

(yj(de) > Also, seven year WMAP observations predict 
that —1.55 < Wde < —0.7 (see, Jarosik et al. [7]). Therefore, the derived model predicts the variation of 
w ( de ) consistent with the seven year WMAP observations (see, Fig. 2). 
Adding ([36]) and ([3?]) . we get the overall density parameter 

Q = Q( rn ) + n(*0 = i _ ^{nDt + d)^^ 1 =l-^A. (38) 

This shows that the overall density parameter depends on the anisotropy parameter A. Since A 
asymptotically tends to provided re < 3 at late times, it follows that O ~ 1 for sufficiently large values 
of t. Thus, the derived model predicts a flat and isotropic Universe for large times, which is consistent 
with the observations. 

Fig. 3 demonstrates the behavior of density parameters in the evolution of Bianchi-I Universe for 
selected values of constants. We observe that initially the ordinary matter density dominates the Universe. 
But later on, the DE dominates the evolution, which is probably responsible for the accelerated expansion 
of the present-day Universe. The overall density parameter curve becomes straight after the dominance 
of DE. Thus, DE may be a possible cause of flatness of the Universe. 

3.2 DE Cosmology for n = 

In this case, integration of ([2]) yields 

a(t) = c 2 e Dt , (39) 

where C2 is a positive constant of integration. 
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0.2 0.4 0.6 0.8 



Figure 3: Density parameters vs t with n = 0.05, D = 2.4, c = 12, b = c\ = 1, u; (m) = 0. 



The metric functions, therefore, read as 



A(t) = aic 2 exp ^Dt - e 

= a 2 c 2 exp |^Z)t - e 

/ 63 
C(i) = a 3 c 2 exp |^£>f - 



1 „-3Dt 



2 -3£>t 



-SDt 



Hence, the model ([TJ becomes 

ds 2 = -di 2 + 4e 2Dt (a 2 e hb ^ 

where h 



3Dt dx 2 + a 2 2 e hb ^ 3m dy 2 + 4e hb ^ 



dz' 



(40) 

(41) 
(42) 

(43) 



-2 

3Dc| ' 

The Hubble parameter, energy density of perfect fluid, DE density and EoS parameter of DE, for the 
model (l43j) are obtained as 

H = D, 



Jm) _ -3(l-Ho<™>) -3£>(l+to< m ))t 
pi — coc 2 e > 



p (de) = 3L) 2 _ p QC -G e ~. _ CoC2 - 



-6Dt 



-3(lW ro >) -3£>(lW m >)t 



(«fe) 



(de) 



P 



-3D 2 - /3 c 2 - 6 e- 6Dt 



c c 2 3(1+, " (m)) w; {m) e- 3D(1+w(m)) * 



(44) 
(45) 
(46) 
(47) 



The above solutions satisfy the equation (|T9j) identically, as expected. 
The spatial volume and expansion scalar of the model read as 



V = 4e 



3'iDt 



(48) 



S 



e = 3D. 

The anisotropy parameter and shear scalar of the model are given by 



(49) 



A = 



2/3oc 2 ~ b 
3D 2 6 



-6Dt 



(50) 



2 o -6 -6Dt 

a = P0C2 e 



(51) 



The density parameters of perfect fluid and DE are as follows: 



Q(m) = c o c 2 



3(lW m) ) 



e -3D(l+w^)t 



3D 2 



(52) 




6Dt C C 2 



3(lW m) ) 



e -3D(l+«i( m ')t 



3D 2 



(53) 



The DP is given by 



q = -l. 



(54) 



Recent observations of SN la ([23] -[32]) suggest that the Universe is accelerating in its present state 
of evolution. It is believed that the way Universe is accelerating presently; it will expand at the fastest 
possible rate in future and forever. For n = 0, we get q = —1 ; incidentally this value of DP leads 
to dH/dt = 0, which implies the greatest value of Hubble's parameter and the fastest rate of expansion 
of the Universe. Therefore, the derived model can be utilized to describe the dynamics of the late time 
evolution of the actual Universe. So, in what follows, we emphasize upon the late time behavior of the 
derived model. 

Fig. 4 shows the behavior of the EoS parameter w^ de > of DE. We observe that u?( rfe ) varies within the 
phantom region and never crosses the PDL. For sufficiently large times, w^ de ^ ~ —1. Therefore, late time 
dynamics of the Universe can be described by DE represented by cosmological constant. 

Further, at late times, we have 



This shows that the ordinary matter density becomes negligible whereas the accelerated expansion 
of the Universe continues with non-zero and constant DE density at late times, as predicted by the 
observations. Also, the Universe remains isotropic at late times. Fig. 5 illustrates the behavior of the 
density parameters during the evolution of the Universe in the derived model. It is observed that the DE 
component dominates the Universe at late times and Q ~ 1. 



pirn) _ Q) p (de) _ ^2 



and 
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Figure 4: vs t with D = 2, fa = cq = c 2 = 1, w (m) = 0. 
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Figure 5: Density parameters vs f with D = 1.2, cq = 3, /3q = C2 = 1, = 0. 



4 Kinematics Tests 

We now study the consistency of the models, discussed in the last section, with the observational param- 
eters. We know that as we look back from our position at r = 0, t = to to some object located at r%, 
we deal with some past moment of time t\ < to- Since r\ and t\ are not directly measurable, we 
measure some physical parameters such as redshift, proper distance, luminosity distance, angular diame- 
ter etc. The various cosmic measures sensitively depend on the parameters of the models, and therefore 
the physical properties of the distant objects are also influenced by these parameters [33l[Mj- In what 
follows, we discuss some important cosmic measures such as lookback time, proper distance, luminosity 
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distance, angular diameter distance and event horizon for the models, discussed in the cases n ^ and 
n = , respectively. 



4.1 Lookback Time-Redshift 

The lookback time is the radiation travel time, taken by a photon emitted by a source at an instant t z in 
the past, to reach an observer at the present epoch to (z = 0). For a given redshift z, the scale function 
a(t z ) is related to a(to) by 

(») 

a{t z ) 

Therefore, the lookback time is exactly the time interval At = to — t z , required by the Universe 
to evolve between these two values of the scale factor. An object at z = 1 emitted its light when the 
Universe was half its present scale (a = 0.5ao). How long ago the light was emitted (the lookback time), 
it depends on the dynamics of the Universe. 

Substituting ([22]) into (f55|) . we get 

t z = t (l + z)- n , (56) 

which can be expressed as 

H (t -t z ) = -[l-(l + z)- n ], (n^O) (57) 
n 

where Hq is the the Hubble's constant in km s _1 Mpc -1 , and its value at present is 72 km s -1 Mpc -1 . 
However, note that Hq is dimensionally similar to the reciprocal of time. The reciprocal of Hubble's 
constant is called Hubble time Tjj : Tjj = Hq~ , where Tjj is expressed in s, and Hq in s^ 1 . If 
Hq is expressed in km s _1 Mpc" 1 , and Tjj in gigayears (1 Gyr = 1 billion years = 10 9 years), then 
Tjj = 977.8 Hq 1 . It is important to note that the age of the Universe can be less than or greater than 
the Hubble time, depending on the values of u;( m ) and u/ rfe ) . 

Taking the limit z — > oo in (|57p . the present age of the Universe (the extrapolated time back to the 
bang) is 

to = (58) 
n 

which, for n = 3/2, reduces to the expression of the standard dust model. 
The lookback time, in terms of the DP, may be obtained as 

H (t - t z ) = -L- \l - (1 + z)-( 1+ «)j . (59) 
1 + q L J 

For small z, equation (f59|) gives 

H {t -t z ) = z- (\ + JV + ( 6 °) 
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Also for q = i, equation ([59|) yields 



2 



-3 



H (to-t t ) = -[l-(l + z)-r\, (61) 

which is used to describe the lookback time in Einstein-de Sitter Universe. In the limit z — > oo, we obtain 
At = (2/3)T H . 

For n = , i.e., q = — 1, we find the following expression for the lookback time: 

Ho(t - t z ) = ln(l + z) . (62) 
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Figure 6: Lookback time versus the redshift z . 



For small z, we have 



Ho(to — t z ) 



1 



2 



(63) 




In Fig. 6, we plot the lookback time as a function of the redshift z for Hq = 100 km s~ Mpc , 
and selected values of q. We observe that the lookback time is larger for negative values of q. Thus, in 
an accelerating Universe, the light from a distant object to reach an observer at the present epoch, would 
take larger time in comparison to the decelerating one. It is maximum in case of q = — 1. 

4.2 Proper Distance-Redshift 

If an observer at r = and t = to receives the light emitted by a source at any instant t located at 
r = r\ with redshift z, then the proper distance between the source and observer is given by 




(64) 



where the radial coordinate distance r%(z) is given by 




(65) 



provided 1. 

Therefore, equation (|64p leads to 




(66) 



In terms of q, the proper distance can be written as 




(67) 



For small values of z, we get 




(68) 
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Taking the limit as z — > oo in equation (j67j) . we get the maximum proper distance as 



d{z = oo) 



H, 



-i 



(69) 



Substituting q = \ into (|67p . we get the following proper distance for Einstein-de Sitter Universe: 



d(z) 



Ho 



For n = , we find 



l-(l + z) : 



1 



d(z) = —z 
-Wo 



(70) 
(71) 
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Figure 7: Proper distance d{z) versus redshift z 



In Fig. 7, we show the proper distance as the function of redshift for H$ = 100 km s _1 Mpc -1 , and 
selected values of q. We observe that the proper distance is larger in case of the accelerating models of 
the Universe. The proper distance curves become almost straight as q gets values closer to —1. There 
is a linear relationship between proper distance and redshift for q = — 1. 

4.3 Luminosity Distance-Redshift 

The luminosity distance of a light source is defined as 

d e =(^y = ri (z)a(t )(l + z), (72) 

where L and i are absolute and apparent luminosities, respectively. r\{z) is the radial coordinate 
distance of the object at light emission, and is given by equation (f65|) . In view of (fM|) . the equation (f72|) 
takes the form 

d e = d(z)(l + z). (73) 
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Substituting (f66|) into (f73|) . the luminosity distance-redshift relation reads as 

1 



[n 



,2-nl 



In terms of q, the luminosity distance-redshift relation can be written as 

1 



d ( 



[(! + *)-(! + * 



,1-91 



Expanding (|75p for small redshifts, after some algebra, we find 



(lr 



1 



z+^(l-q)z 2 + 



(74) 



(75) 



(76) 



The luminosity distance for the Einstein-de Sitter Universe (q = |) is given by 

2 



do 



Hn 



;i + z)-(i + z) 



(77) 



For q = 1 , flo^i = z , which shows linear relationship between luminosity distance and redshift. 
For n = , i.e., q = — 1, we find 

d e = -±-(z + z 2 ). (78) 
-Ho 
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Figure 8: Luminosity distance (c^) versus redshift z 



Luminosity distance as a function of the redshift is shown in Fig. 8 for Hq = 100 km s _1 Mpc -1 , 
and selected values of q. We observe that all curves start with the linear Hubble law {H^di = z) for 
small z, but then only the curve for q = 1 shows linearity all the way. We also note that for the same 
redshift, the luminosity distance is larger for smaller values of q. The rate of increase of di is fastest in 
case of q = — 1. All the curves merge for small z. 
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4.4 Angular Diameter Distance-Redshift 

The angular diameter of a light source of proper distance d{z) at r = r% and t = t\ , observed at 
r = and t = tQ, is given by 

«,)(!+,)'_ 
ria(ti) d£ 

The angular diameter distance is defined as the ratio of the source diameter to its angular 

diameter, and therefore 



Inserting (J7U) into (J80]) , we find 

cU = 



(n 



na(ti) = dg(l + z)~ 



In terms of q, the angular diameter distance can be written as 

1 



d A 



qH 



(l + z)- 1 - (1 + z)- 



(1+9) 



For the Einstein-de Sitter Universe, we get 

2 



d A 



H, 



(l + z)" 1 -(l + z) 



(80) 
(81) 

(82) 

(83) 



The extremum value of angular diameter distance occurs at z c = (1 + q)" — 1. Therefore, for the 
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Figure 9: Angular diameter distance (d A ) in units of H 1 versus redshift z. 



Einstein-de Sitter Universe, d A is extremum at z c = |. 



For n = , i.e. 



T, we find 



d A 



1 

77,, V i + •- 



(84) 



In Fig. 9, we show the angular diameter distance versus redshift for selected values of q in the units 
of Hq 1 . We observe that each of the angular diameter distance curve has one point of maxima. For 
all models, the angular diameter distance increases with increasing z , reaches its maximum value at a 
given z c , and eventually begins to decrease. Usually d A has a maximum (or minimum) for some z = z c . 



15 



4.5 Event Horizon 

The event horizon r_e is the limit of the proper distance such that a light signal sent by a distant source 
at r = n, t = ti , will never reach an observer at r = 0, t = to if n > Te- In other words, the two 
observers having a proper distance greater than te , can never communicate each other. Mathematically, 
the event horizon is given by 

, . f°° dt 

r E = a(t ) / -777 • (85) 



to 



a(t) 



Substituting (J32J) into ([85]), we get 



For n = , i.e., q = — 1, we find 



= tft • (86) 

(1 - n)H 



r E = ^~- (87) 



We observe that an event horizon exists in the models where n < 1 or q < . It follows that event 
horizon does not exist in the Einstein-de Sitter Universe. Equation (|87p gives the event horizon, where 
no observer beyond a proper distance re at t = to can communicate with another observer. 



5 Concluding Remarks 

We have presented singular (n ^ 0) and non-singular (n = 0) Bianchi type-I models of Universe by taking 
into account dynamical DE along with the usual matter fluid content. The models describe shearing, non- 
rotating and expanding Universe. The anisotropy of the Universe goes off during the cosmic evolution. 
The analysis of the solutions reveals that the Universe achieves isotropy and flatness after the dominance 
of DE. Thus, one may hope that DE would be a plausible explanation to the unsolved problems of 
cosmology. The values of EoS parameter of DE have been found in the range predicted by seven year 
WMAP observations. Thus the derived models are consistent with the WMAP observations and hence 
are physically viable. 

The explicit expressions for the observational parameters such as the lookback time, proper distance, 
luminosity distance, angular diameter distance and event horizon have been obtained in the power-law and 
exponential-law cosmologies. The graphical analysis reveals that the lookback time, proper distance, lu- 
minosity distance and angular diameter distance have larger values in accelerating models of the Universe. 
The maximum values have been observed in the case q = —1. The event horizon exists in the accelerating 
models of the Universe. It may be noted that the analytical expressions of the observational parame- 
ters in power-law and exponential-law models have been discussed in many works \35 \ 1361 [37] . However, 
in addition the present work presents the graphical analysis of the parameters, and demonstrates their 
behavior in accelerating models in contrast with the decelerating models. 

It is worthwhile to mention that although any constant DP model is not capable of representing the 
whole or a long period of history of the Universe at one go, it can still be considered as an approximation 
for representing a particular period of the history of the Universe, for instance, the vicinity of the present 
time of the Universe |23[ l32j. Thus the models presented in this work can be utilized for this purpose. For 
instance, the non-singular model (exponential-law model) successfully describes the late time evolution 
of accelerating Universe with DE and puts forward cosmological constant as a candidate of DE. 
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